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The low temperature scanning tunneling microscopy spectra in the underdoped regime is ana- 
lyzed from the perspective of coexisting d-density wave (DDW) and d-wave superconducting states 
(DSC). The calculations are carried out in the presence of a low concentration of unitary impurities 
and within the framework of the fully self-consistent Bogoliubov-de Gennes theory, which allows 
local modulations of the magnitude of the order parameters in response to the impurities. Our 
theory captures the essential aspects of the experiments in the underdoped BSCCO at very low 
temperatures. 
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A fundamental tension in the field of high tempera- 
ture superconductors is the notion of a competing order 
parameter in the underdoped regime, which can provide 
a natural explanation of why the superconducting dome 
exists and shed light on the nature of the pseudogap. 
While a charge ordered state is a candidate [l] , one of us 
has proposed, and extensively studied, a new order pa- 
rameter, which results in circulating currents arranged 
in a staggered pattern(DDW) Q . Many experiments are 
consistent with this order parameter, as demonstrated 
in studies of the superfluid density, the polarized neu- 
tron scattering, the Hall number in pulsed 60 T mag- 
netic field, the angle resolved photoemission spectroscopy 
(ARPES), the lack of specific heat anomaly at the pseu- 
dogap temperature, the transition temperature in multi- 
layer cuprates, and the infrared Hall angle measurements 
|sj|. So far the clinching direct experiment, the polar- 
ized neutron scattering, has remained suggestive |^ be- 
cause of the difficulty of detecting weak signals from the 
small magnetic moments generated by the circulating or- 
bital currents. Recently, another novel experimental test 
has been suggested that takes advantage of the spin-orbit 
coupling in the DDW state 0- 

Here we turn to the intriguing scanning tunn eling mi- 
croscopy (STM) measurements [& 0, H, M, OH ED- In 
spite of numerous theoretical analyses of this problem 
El El El El El El El El El El, no compre- 
hensive theoretical picture has yet emerged, although cer- 
tain aspects are captured by some of them. For example, 
earlier measurements in the slightly overdoped samples 
at very low temperatures have been elegantly explained 
in terms of a quasiparticle scattering interference model, 
named the octet model HQ]. At the same time, an in- 
terpretation in terms of dynamic charge fluctuations has 
also been advanced The focus here, however, is 

on an extensive set of experiments as a function of dop- 
ing in BSCCO at very low temperatures 0. The exciting 
finding of these experiments is the emergence of a new 
order, present along with the d-wave superconductivity 



(DSC). The salient signature is a sudden development of 
a relatively non-dispersive incommensurate wave vector, 
q*, at higher energies in the underdoped regime. 

In this Letter we explain the experiments by adopting 
a view that is orthogonal to the notion of charge order 
and consider DDW. At first sight, this would seem im- 
possible, as STM is not sensitive to currents. The key 
to this puzzle lies in the presence of disorder (inevitable 
in these materials) that leads to variations of the charge 
density, which in turn can scatter the quasiparticles of the 
DDW, revealing indirectly the DDW order in the form 
of an interference pattern. The idea is not as surpris- 
ing as it may seem, because, after all, this was precisely 
the basis of the successful octet model. Clearly, DSC 
corresponds to off-diagonal long range order rather than 
diagonal long range order and therefore cannot directly 
affect STM. The general message is that a large class 
of ordered states can leave their characteristic signatures 
via the scattering interference of their distinctive quasi- 
particle spectra. Our analysis is broadly consistent with 
the conclusions drawn from the recent comparisons of the 
autocorrelation function of the ARPES spectra and the 
STM spectra HHH. 

A single-impurity T-matrix calculation in the coexist- 
ing DDW and DSC state El 0] was unable to recover 
the salient feature of the experiments, while pin-pointing 
certain important aspects . We shall show that the fault 
lies with the method and not intrinsically with the idea 
that a new order arises in the low temperature under- 
doped regime, namely the coexisting DDW and DSC or- 
der. The problem is that the T-matrix approach neither 
allows the amplitude of the order parameter of the DSC 
to modulate in response to the impurity potential nor al- 
lows a proper treatment of the current conservation in the 
DDW state in the presence of impurities [H E|| . More- 
over, it excludes not only spatial structures but also the 
interference of the impurities. 

To understand the effect of various competing orders 
on the local density of states (LDOS), all we need to know 
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are the corresponding wave functions in the presence of 
finite concentration of impurities. We accomplish this by 
choosing the simplest Hamiltonian that yields both DSC 
and DDW, and coexisting DDW and DSC, states within 
the Bogoliubov-de Gennes mean field theory. We use 
this Hamiltonian solely as a crutch to generate the states 
with the broken symmetries that we wish to study. In 
general, deep within the superconducting dome, that is, 
at very low temperatures, there is good evidence of the 
existence of quasiparticles. The experiments that we ad- 
dress here are precisely in the regime in which our theory 
is expected to be valid. The corollary to these observa- 
tions is that a straightforward finite temperature exten- 
sion of our theory cannot be applied to the experiments 
of Ref. because the simple picture of quasiparticles 
may not be valid at higher temperatures and above the 
superconducting dome. 

The Hamiltonian is TL = K + H ul t + Hdis, where K = 
~J2ij, a ( t v c ia c ja + h.c). The operator c\ a (c ia ) creates 
(destroys) an electron at the site i with spin a, and Uj 
is the hopping matrix element to the nearest (t) or to 
the next-nearest (— 0.2732: ) neighbor J^. The interacting 
part consists of 
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which is identified by two parameters J and W |28| (S^ 
is the spin operator at the site i and rii a is the den- 
sity operator). Disorder is introduced through Hdis = 
Y^ia (V* ~ ^) n ia, Vi, at each site i, is an independent 
random variable, which is either +Vq (repulsive) with a 
probability n lmp (impurity concentration) or zero, and fi 
is the chemical potential of the system. 

We solve the above model within the self-consistent 
framework of BdG mean field theory, which amounts to 
decoupling the interaction terms, resulting in a TL c s that 
leads to the local order parameters: DSC order (Aj), 
DDW order (xi) an d also the local Hartree and Fock 
shifts. The details can be found in Ref. [2(| . The chem- 
ical potential /i is obtained by fixing the average density 
of electrons, (n) = J2i(ni)/N, at the desired value (N 
being the system size). 

We concentrate on the results at T = for parame- 
ters J = 1.6, W = 0.6 in units of t (all energies will be 
expressed in this unit and all lengths will be expressed 
in terns of the lattice spacing), and (n) = 0.9 (equiva- 
lently 10% hole doping) on square lattices with a unit 
cell of size N = 30 x 30. Larger unit cells result in neg- 
ligible improvements at the cost of computer time; the 
dimensionality of the matrices for which all the eigenval- 
ues and the eigenvectors must be repeatedly obtained is 
2N x 2N, that is, 1800 x 1800. These parameters result 
in xo = 0.33 and A = 0.12 in a disorder free system in 
the coexisting DDW and DSC phases, where \o is de- 
fined through Xk = Xo(cosfcic — cosfc y )/2 and similarly 



for Ao. Results for systems with solely DDW or DSC 
order will also be discussed for slightly different param- 
eters: J = 1.4, W = 0.5. We choose V = 100, close to 
the unitary limit, and results are averaged over 8-10 dif- 
ferent realizations (8 for the DDW and the DSC orders 
by themselves, but 10 for the coexisting DDW and DSC) 
of the random potential. For better statistics we use the 
repeated zone scheme jl^, in which a super-cell is 
constructed by replicating the 30 x 30 unit cell 10 times 
along the x-direction and 10 times along the y-direction. 
This provides a denser set of energy eigenvalues. 



(Jt,7 t) 



(a) 



(7t,7 l) 



(b) 



(-Ji,-Ji) m = _ .3i5 Hi-ji) co = —0.245 

(71,71) (7t,7l) 



(O 



(d) 



(-Tt-Ti) = +Q.315 (-Ji-Ti) co = +0.280 

(Tt.Tl) 



(e) 



(71,7 1) 



(0 



(-7E-7C) m _ +0 .070 



(-JI-JI) m _ aoo(J 



FIG. 1: Density plots on a linear scale of the FT-LDOS 
within the first Brillouin zone in the k x — k y plane. The 
panels (a) through (f) correspond to different energies, uj, for 
the coexisting DDW+DSC phase. Here m mp — 0.01 and Vo = 
100. Peaks at q* ~ 0.15(27r) are clear in panels (a) through 
(d). At smaller energies, panels (e) and (f), g*-peaks are 
absent and a broad peak appears near q = 0, whose width 
decreases with increasing The strength of the peak at 
(0, 0) is reduced for the clarity. 

In the presence of impurities, both orders become spa- 
tially inhomogencous. However, at n- lmp = 0.01, the (dis- 
order) averaged DDW order degrades much less (xi ~ 
0.95xo) compared to the DSC order (Aj w 0.55A ), con- 
sistent with earlier findings |2(|. We define the Fourier 
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transform iV(q, uj) = J2 T e ,q r iV(r, of). 
the LDOS at a site r with energy u>. 



where N(r, uj) is 
The disorder av- 



eraged power spectrum P(q, uj) = \N(q,uj)\ 2 /N is calcu- 
lated and then yj -P(q, uj) is compared with experiments, 
as this is the theoretical measure |l3j of the Fourier trans- 
form of the experimentally measured dI(V, r)/dV, where 
/ is the tunneling current, and V is the bias voltage 
{uj = eV). 

The results in the coexisting DDW+DSC state are 
shown in in Fig. ^ The intensity of the peaks along 
(0, ±1) and (±1, 0) appears as a generic and robust prop- 
erty, more precisely at q* ~ (0.15)27r for the chosen set of 
parameters. The q* peaks occur for \uj\ > Ao. At lower 
energies the FT-LDOS profile looses periodic q* modula- 
tions and a rather broad peak occurs in its Fourier trans- 
form at q w OfSee Figs. 1(e) and 1(f)). Unitary impurity 
resonances |81j dominate this regime (near uj ss 0) and 
wash out the globally periodic LDOS modulations aris- 
ing from interference. Unitary scatterers do not affect 
the spectrum at energies larger than Ao, and hence the 
LDOS modulations are preserved. 

The finite system, 30 x 30, allows a q-resolution of 
0.033, which is rather coarse compared to the real data 
(the super-cell does not introduce any additional inde- 
pendent wave vectors due to the periodicity) . Along with 
the strong peak at q = mentioned above, the resolu- 
tion problem prevents us from observing the low energy 
modulation q\ , in the terminology of Ref. 9] ; it is hidden 
by the broad peak around q = 0, whose width decreases 
with \uj\. On the other hand, q^ is quite visible (at least 
for lu > 0). It is known from Ref. |9j that the intensities 
of (72, <?3, ?6i 97 relative to q\ and or f all with decreasing 
doping and q^ is not seen at all [3, [^J. This is con- 
sistent with our calculation, which was performed in the 
underdoped regime. 

Because DSC order is already weak in the pure system 
for our chosen parameters, and it rapidly gets weaker 
with impurities, it is important to study profiles similar 
to Fig. 1 in the DDW phase without coexisting DSC order 
with 7ii mp = 0.01. The results are presented in Fig. 2. 
The features are similar to those of Fig.l, but the q* 
peaks first emerge for a little smaller value of \oj\ (see, 
Fig. 13) • The spectra shown in Fig. [21 are generic, and 



similar results are obtained for m, 



0.02. We have 



also studied our model in the DSC-only state in the same 
underdoped region and have not found q* peaks. 

An intuitive explanation of q* is as follows. It is due to 
scattering between the tips of the hole pockets and the 
sharp rise of the pseudogap. In the coexisting DDW and 
DSC state, the low-energy constant energy contours are 
"bananas" , essentially from the pure DSC order, except 
that they are doubled up because of the broken trans- 
lational symmetry However, the disorder considerably 
smears out these scattering events. As the energy in- 
creases these bananas coalesce and, then on, the scatter- 
ing is between the tips of the hole pockets of DDW. Since 



DDW is less affected by potential scattering than DSC, 
the signature becomes robust. 
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FIG. 2: FT-LDOS as in Fig.l but for DDW (%o « 0.28). 
DSC order is forced to zero on all sites. Qualitative features 
are similar to Fig.l for a wide range of u. 
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FIG. 3: Top panel: dispersion of q* (in units of 2n) with u) in 
an impure DDW+DSC phase. Bottom panel: the same as in 
top panel but for the DDW phase. The error bars are about 
the size of the symbols. 

The dispersion of the g*-peaks is shown in Fig. 01 The 
dispersion in the top panel of Fig.[3|is somewhat stronger 
than in Ref. but gets weaker in the DDW-only phase 
(bottom panel) . Similar results are also found for ni mp = 
0.02. The value of q* ~ (0.15)2tt is due to the chosen set 
of band parameters, which will change with a different 
choice. 

We have repeated our calculations with impurity 
strengths Vq = 10 and Vq = 1. The results are essen- 
tially the same for Vq = 10, while for Vq = 1 (non- unitary 
scatterers) the spectra for negative energy are somewhat 
different, though the key features are preserved. Thus 
the robustness of the results found here is heartening. 
We emphasize, however, that the effects of inhomogene- 
ity and impurity interactions can be very subtle 30 . 

The disorder averaged DOS, N(u), is shown in Fig. 0] 
We see that N(uj) is asymmetric in the DDW phase, 
both with and without DSC order. The asymmetry in 
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FIG. 4: Impurity averaged normalized DOS, N(uj), for co- 
existing DDW+DSC for ni mp = 0.01. The spectrum is asym- 
metric and low energy humps appear around Ao. The inset 
corresponds to DDW (xo ~ 0.28) and DSC (A ~ 0.196), 
with n imp = 0.01, J = 1.4, and W = 0.5. 



the spectrum has the same origin as the asymmetries in 
Figs. Hand 12 In the coexisting phase, the remnants of 
the superconducting density of states peaks can be seen, 
below which modulation of LDOS dies out. The spec- 
trum is tantalizingly similar to the experimental results 
of Ref. 0. 

To summarize, the inhomogeneous phase of coexist- 
ing DDW and DSC captures the essential experimental 
features seen in the underdoped regime at low tempera- 
tures, namely the emergence and the nature of the rel- 
atively nondispersive q* peaks seen only at higher en- 
ergies; in contrast, the DSC order alone is inadequate. 
From our perspective the emergence of q* is not due to 
an explicit charge order, fluctuating or otherwise, but 
due to the current modulations of the DDW. It would 
be interesting to see if q* extends to lower energies with 
purer samples. The interpretation of the experiment 
in Ca2- x Na x Cu02Cl2 is an open issue. 

This work was supported by the grant NSF-DMR- 
0411931. We would like to thank C. Bena, J. C. Davis, 
J. Hoffman, H. -Y. Kee, K. McElroy, S. Kivelson, and 
D. J. Scalapino for helpful comments. We acknowledge 
the computational support of H. U. Baranger's group at 
the Duke University, where part of the calculations were 
carried out. 



[1] S. A. Kivelson et al. Rev. Mod. Phys. 75, 1201 (2003). 
[2] S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. 

Nayak, Phys. Rev. B 63, 094503 (2001). 
[3] S. Tewari, H. -Y. Kee, C. Nayak, and S. Chakravarty, 

Phys. Rev. B 64, 224516 (2001); S. Chakravarty, H.-Y. 

Kee, and C. Nayak, Int. J. Mod. Phys. B 15, 2901 (2001); 



S. Chakravarty, C. Nayak, S. Tewari, and X. Yang, Phys. 
Rev. Lett. 89, 277003 (2002); K. Hamada and D. Yosh- 
ioka, Phys. Rev. B 67, 184503 (2003); H. -Y. Kee and Y. 
B. Kim, Phys. Rev. B 66, 012505 (2002); S. Chakravarty, 
Phys. Rev. B 66, 224505 (2002); H. -Y. Kee and Y. B. 
Kim, Phys. Rev. B 66, 052504 (2002); S. Chakravarty, 
H. -Y. Kee and K. Voelker, Nature 428, 53 (2004); S. 
Tewari et al., Phys. Rev. B 70, 014514 (2004). 
[4] H. A. Mook et al. Phys. Rev. B 69, 134509 (2004); H. A. 

Mook et al. Phys. Rev. B 66, 1 44513 (2002) 
[5] C. Wu, J. Zaanen, S.-C. Zhang, eond ma! .05055 11 
[6] J. E. Hoffman et al. Science, 297, 1148 (2002). 
[7] C. Howald et al. 14533 (2003). 
[8] K. McElroy et al. Nature, 422, 592 (2003). 
[9] K. McElroy et al. Phys. Rev. Lett. 94. 197005 (2005) and 
|cond-m at /0404005] . 
[10] M. Vershinin e t all Science, 3 03, 1995 (2004). 
[11] A. Fang et al, cond-mat/0404452 A. Kapitulnik et al, 

cond-mat/0407743 
[12] J.M. Byers, M.E. Flatte and D.J. Scalapino, Phys. Rev. 
Lett. 71, 3363 (1993); M.E. Flatte, and J.M. Byers, Phys. 
Rev. Lett. 80, 4546 (1998). 
[13] L. Capriotti, D. J. Scalapino, and R. D. Sedgewick, Phys. 

Rev. B 68, 014508 (2003). 
[14] Q.-H. Wang and D.-H. Lee, Phys. Rev. B 67, 20511(R) 
(2003). 

[15] C. Bena, S. Chakravarty, J.-P. Hu and C Nayak, Phys. 

Rev. B 69, 134517 (2004). 
[16] C. Bena, S. Chakravarty, J.-P. Hu and C Nayak, 

cond-mat/0405468 
[17] D. Podolski, E. Demler, K. Damle, and B. I. Halperin, 

Phys. Rev. B 67, 094514 (2003). 
[18] C. -T. Chen, and N. -C. Yeh, Phys. Rev. B 68, 220505 

(2003) . 

[19] A. Polkovnikov, S. Sachdev, and M. 

388-389 19 (2003). 
[20] H. -D. Chen, O. Vafek, A. Yazdani, S. 

Rev. Lett. 93, 187002 (2004) 

[21] H. C. Fu, J. C. Davis, and D. -H. Lee, |cond-m'a t/0403001 
[22] T. Pereg-Barnea and M. Franz, Phys. Rev. B 68, 

180506(R) (2003). 
[23] U. Chatterjee et al. cond-mat/0505296 The definition of 

q* in this paper is different from ours. Our q* exists only 

at high energies, as in Ref. |9j. 

[24] K. McElroy et al. cond-mat/0505333 

[25] A. Ghosal, M. Randeria, and N. Trivedi, Phys. Rev. B 

63, 20505(R) (2000). 
[26] A. Ghosal and H.-Y. Kee, Phys. Rev. B 69, 224513 

(2004) . 

[27] Highly complex BdG equations restricted us to first three 
non-trivial band structure parameters from M. R. Nor- 
man, M. Randeria, H. Ding, J. C. Campuzano, Phys. 
Rev. B 52, 615 (1994). The precise value of q* , but not 
its existence, depends sensitively on the band parameters 
16]. 

[28] The J and W terms contribute to both DDW and DSC 
at the mean field level. We keep both to tune them inde- 
pendently. See Ref. |2^] as well as C. Wu and W. V. Liu, 
Phys. Rev. B 66, 20511(R) (2002) for further details. 

[29] A. Ghosal, C. Kallin and A. J. Berlinsky, Phys. Rev. B 
66, 214502 (2002). 

[30] L. Zhu, W. A. Atkinson and P. J. Hirschfeld, Phys. Rev. 
B 69, 60503(R) (2004). 

[31] S. H. Pan et al. , Nature, 403, 746 (2000). 



Vojta, Physica C 
-C. Zhang, Phys. 



T.Hanaguri et al. Nature 430,1001 (2004). 



